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We consider conductance mapping of systems based on the two-dimensional electron gas with 
scanning gate microscopy using two and more tips of the atomic force microscope. The paper 
contains results of numerical simulations for a model tip potential with a proposal of a few procedures 
for extraction and manipulation the ballistic transport properties. In particular, we demonstrate 
that the multi-tip techniques can be used for readout of the Fermi wavelength, detection of potential 
defects, filtering specific transverse modes, tuning the system into resonant conditions under which 
a stable map of a local density of states can be extracted from conductance maps using a third tip. 


I. INTRODUCTION 

Conductance ( G ) of open systems with the two- 
dimensional electron gas (2DEG) in semiconductor 
heterostructures at low temperatures is determined by 
scattering of the Fermi level electrons. Using the scan¬ 
ning gate microscopy (SGM) technique 1 one probes 
the properties of the electron transport in mesoscopic 
devices by variation of the potential landscape for 
the Fermi level electrons with the charged tip of the 
atomic force microscope (AFM) moving above the sur¬ 
face of the sample. The potential of the charge at 
the tip is screened by the two-dimensional electron 
ga^ ? which produces a short-range form of the effec¬ 
tive tip potential for the Fermi level electrons. The 
SGM technique was widely used in the studies of the 
ballistic electron transport, in particular to visual¬ 
ize the electron trajectories as deflected by external 
ma gnet ic fielcP, electron branching of the Coulomb 
flowS^l including an evidence of the quantum Braess 
paradosP, formation of a quantum ring potential with 
a controlled number of modes in each arnP, investi¬ 
gation of el ectron backscattering with quantum point 
contactiPP formation of Coulomb islandtP and map¬ 
ping the local density of states (LDOS^™. Besides 
the studies of the ballistic flow, the tip potential is also 
used for the studies of the charge flow in conditions of 
the Coulomb blockade in quantum dots 1 3. The role of 
the tip is then to tune the chemical potential of the 
confined electron system into the transport window 
defined by the Fermi energies of the source and drain 
electrode^ 1 . 

In this paper we consider possible applications of 
the SGM technique for the studies of the ballistic 


transport using two or more tips instead of a single 
one. We show that the double tip system can be used 
to: i) measuring the Fermi wavelength, ii) mapping 
the potential defects in the channel iii) mode filtering, 
iv) detection of localized resonances in an experimen¬ 
tal implementation of a stabiliation method, v) tuning 
the system to resonant conditions when LDOS can be 
read-out with a third tip. 

Usage of several probes was implemented a few 
years ago for the scanning tunneling microscopy 
(STM) studies of the sample surface. A version of 
the STM technique using several independent tips 
was used in particular to perform the four point 
measurements-^ for determination of the surface prop¬ 
erties. 


II. MODEL 

We consider the ballistic transport at the Fermi 
level electrons in systems based on a two dimensional 
electron gas (2DEG). We neglect the electron-electron 
interactions and we use the effective mass Hamiltonian 
of form 

{ 2m”-fr ^ Vext( x i 2/)} v) = E Fip{x,y), (1) 

where m e s = 0.067mo is the effective mass of GaAs, 
E-p is the Fermi level energy, and T4 x t contains all 
external sources of the electrostatic potential (eg. po¬ 
tential of the tips V)ip)- We assume that the potential 
of the tip is given by short-range Lorentzian potential 
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This form of the potential was evaluated in previous^!, 
self-consistent Schrodinger-Poisson calculation. In the 
following we assume that all the tip potentials have 
the same width and amplitude and that the distance 
between them can be changed. The shape of the dis¬ 
cussed devices - tailored from the sample containing a 
two-dimensional electron gas is described by hard-wall 
boundary conditions. 

In order to solve the scattering problem we use a 
finite difference implementation of the Transparent 
Boundary Method (TBM^PMI 

For the boundary 
conditions in the input lead we use the standard ap¬ 
proach with the wave function given by superposition 
of incoming and outgoing (reflected) transverse modes 


-Minput 

V’input (x,y) = [ake ikx Xk PU \y) ( 3 ) 
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where Mi nput is the number of current propagating 
transverse modes Xfc nput the i n P u t lead, a*, and rx 
are the incoming and reflection amplitudes. For the 
output lead we assume that the wave function is given 
by formula 
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^output Or, y) = ^ kx xT PUt (y), (4) 

k =1 

where M outpu t is the number of transverse modes 
^output ^-j le ou tp U t lead and tk is the outgoing ampli¬ 
tude. The transverse modes were calculated with the 
method presented in RefP^. Matching the boundary 
conditions © and 0 with the wave function calcu¬ 
lated inside the device one finds the solution for the 
scattering problem. After solution of the Eq. Q (for 
details see 1 ^*) for each incoming mode i the conduc¬ 
tance of the system is calculated from the transmission 
probability T) using the zero-temperature Landauer 
formalism 


-^input 
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with Go = ^f-. 


In all the cases considered below we assume that 
both leads have the same width, thus M = Mi nput = 
-^output- We choose the discretization grid Ax = 
Ay = 4nm. In further discussion we will refer to the 
local density of states which is defined as the sum 
of electron densities incoming from the left and right 
leacP. 


III. IMPURITIES MAPPING 


Let us start our discussion with the device presented 
in Fig. [l] We consider a long (lOOOnm) and narrow 
(80nm) channel. For the Fermi energy E-p = 2.5meV 
there is only one current-propagating transverse mode 
in the channel (M = 1). We consider a double tip sys¬ 
tem assuming that tips move along the x axis. The 
center of the system of the two tips is denoted by x t i p 
(see Fig. [I]). Both tips are separated by a distance 
of dtip- We assume Ut\ p = 5meV and wti p = 15nm, 
For a single tip with these parameters the conduc¬ 
tance of the system is reduced to G = 0.12Go- In Fig. 
[2ja) we show the result for the conductance G of the 
system as a function of Xu p and the intertip distance 
G?ti p for a clean channel. The d t ; p dependence of con¬ 
ductance reveals a series of resonances separated by 
distance A^/2 ~ 60nm, which are related to forma¬ 
tion of standing waves between the tips. Thus with 
the two tips one can probe the Fermi wavelength. 


^tip scan line 



L = lOOOnm 


Figure 1: Sketch of the system considered in Sections 
uni and nvi We consider a channel of width W — 80 
nm and the system of two tips with the same amplitude 
[/tip = 5meV and width Wti p = 15nm The length of the 
computational channel is L = lOOOnm. The tips are above 
the axis of the channel (red dashed line) and are separated 
by a distance d t i P - The center of the system of two tips is 
located at position Xti P - 

Let us now consider the system with a single po¬ 
tential defect present in the middle of channel. The 
defect is modeled with Eq. ([2]) with Z7i mp = 2meV 
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and wimp = lOnm (see the potential profile below the 
Fig- §b)). Now the SGM image reveals new features 
- the resonance lines bend and oscillate in the vicinity 
of the impurity. Far from the impurity the resonance 
lines return to the same position as in Fig. [2](a). Note 
that the first resonance line (around the d t ip ~ 90nm) 
resemble the potential profile inside the channel. In 
Fig. ,d) we show that the first resonance line fol¬ 
lows the potential profile of the defect. This effect can 
be explained in terms of the semi-classical WKB ap¬ 
proximation. For slowly varying potential U(x) and 
E-p > U ( x ) the wave vector of propagating wave func¬ 
tion e ±lk ^ x is given by k F = ^/2 m e s (E F — U(x))/h. 
Now if U(x) > 0 the phase of the wave propagat¬ 
ing through the potential hill will be delayed by some 
value Aip in comparison to the phase of the unper¬ 
turbed system. Thus in order to restore the resonance 
for standing wave one must increase the distance be¬ 
tween the tips, by the value which compensates for 
the delayed Aip phase. For U(x) < 0 the distance 
between the DT has to be decreased (see Fig. id)). 

The system of a single tip and a repulsive defect is 
equivalent to the double-tip system, hence the position 
of the defect can be resolved by a scan of a single-tip 
by a shift of the energy lines. Nevertheless, in order 
to extract the profile of the potential of a defect one 
needs two tips and not a single one. 

In Fig. [3] we show the results for the channel with 
two impurities (same as in Fig. |2'c)) but for higher 
energies E F . At higher energies (Fig. la)] the res¬ 
olution of the images is reduced with the resonance 
lines getting closer and wider, but the image still al¬ 
lows one to map the defect potential distribution along 
the channel. Note that in Fig. |3](a) we notice the 
resonances of Fig. [2jc) for both the lowest subband 
(resonance sequence starting near du p — 30 nm) and 
the second transverse mode (the sequence starts for 
d t ip — 100 nm). For E F = 7.5meV (see Fig.|3j(b)) 
the image looses the features visible in Figs. [2pbut 
the W shaped line clearly indicates the position of the 
defects. 

For completeness we considered a clean channel 
with the width of the tip potential increased to i«tip = 
50nm. In Fig. [4] we plotted the conductance as a 
function of the intertip distance du p and the height 
of the tip potential. For Uu p = 2 meV the system 
is transparent for the electron flow, formation of the 
resonances as discussed above for wu p = 15 nm ap¬ 
pear higher in the energy and for U t i P >3.5 meV the 
ballistic electron flow is blocked. For larger values of 
Uti p formation of a quantum dot supported by the tips 



Figure 2: Results for Ep = 2.5meV (one band transport) 
a) The DT scan obtained for a clean channel, b) Added 
impurity in the middle of channel with amplitude U lrnp = 
2meV and of width tUi mp = lOnm modeled with Eq. |2j. 
c) the same but for two impurities separated by distance 
160nm and d) obtained for two impurities of different sign 
±2meV separated by distance lOOnm. Plots below show 
the potential profile in the channel for each case (b-d). 
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Figure 3: The same as in Fig. [2jc) but for higher Fermi en¬ 
ergies. a) Results for Ep = 4.5meV (two band transport) 
and b) Ep = 7.5meV. 


should be expected with the transport dominated by 
the Coulomb blockade, which is however, outside the 
scope of the present work. The best resolution of the 
scan as a function of da p is obtained just below the 
cutoff of the ballistic transport. 

To conclude this Section, we find that double tip 
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Figure 4: Conductance of the system as a function of the 
tip amplitude C/ti P and distance d t i P between the tips ob¬ 
tained for Ep = 2.5meV and width of the tip potential 
ictip = 50nm. 


system can be used to read out the Fermi wavelength 
and the potential profile along the channel. 
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Figure 5: a) The conductance G of the long channel of 
width W = lOOnm as a function of d t i P and Ep. Dashed 
lines deposited to the (a) shows the energy spectrum a two- 
dimensional quantum well of length dti P — 2n;ti P . The shift 
by 2u;ti p accounts for the finite size of the tip potential, 
(b-f) Electron density \ip\ 2 obtained for values of (dti p ,Ep) 
pointed by arrows. 


IV. TRANSVERSE MODE FILTER 

In Fig. [3^a) we observed the resonances of the first 
and second subband (the red lines near d t i P = 100, 
200 and 300 nm ). These resonances occur separately 
or in overlap with the lines of the lowest resonance de¬ 
pending on the values of dtip and .T t i p . Thus choosing 
a proper value of d t i P and cc t i p one can make a device 
which will filter out specific modes. For the proof of 
principle, we will restrict our consideration to the first 
and second mode. The resonance lines are plotted in 
Fig. [5ja) as functions of the distance between the tips 
and the Fermi energy for the clean channel (as in Fig. 
|T|) but with the width increased to W = lOOnm. At 
low Fermi energies only the first transverse mode is 
present in the conductance, and the double tip stops 
the transport unless the Fermi energy coincides with 
the resonances localized between the tips. For higher 
values of Ep the resonance lines of the second mode 
appear and at lines in the (dup^F) plane which in¬ 
tersects with the first mode lines. The images (b-f) 
in Fig. [s] show the electron density for work points 
marked by arrows in Fig. [Sja). We can see that for 
a given energy Ep one may filter out a specific trans¬ 
verse mode, by changing the dti p distance. 


In Fig. [5ja) we plotted the energy spectrum of 
an infinite quantum well of width 100 nm and length 
d t ip ~ 2 w t ip. The conductance resonances follow ex¬ 
actly the energy spectrum for a infinite quantum well. 


V. DETECTION OF LOCALIZED 
RESONANCES 

Let us now consider the system depicted in Fig. 
[0] with a channel of width lOOnm connected to the 
quantum cavity of length 500nm and width W. In 
order to get the best resolved images we set the Ep = 
2meV (single subband transport within the channel). 
In Fig. |7ja) we show the conductance of the system as 
a function of the dti p and width W of the cavity. We 
kept the center of the double tip fixed in the middle 
of the resonant cavity. For small values of W around 
lOOnm the resonance lines are very similar to those 
of Fig. [BJa). For a changed width of the cavity the 
cavity-localized energy levels vary as oc W~ 2 . Thus 
changing the width W we should expect behavior of 
the resonance lines similar to the ones of Fig. [sja) . 
Large values of W lead to a complex behavior of the 
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resonance lines but with well distinguishable patterns 
of X-shaped lines which appear with a period of Xf 
along dtip axis. 



500nm 


Figure 6: Sketch of the second system considered in this 
paper. Long channel of width lOOnm is coupled to the 
resonant cavity of variable width W and length 500nm. 
The same as before we have two tips in the system one in 
the left lead and second in the right lead. Both tips have 
the same parameters of ?7 tip = 5meV and Wti p = 15nm 
and are separated by distance d t i p . Red line corresponds 
to the third tip scan discussed in the last paragraph. 

In Fig. [7^b) we show the conductance of the system 
without the tips. The conductance contains a series 
of sharp resonances of a Fano type corresponding to 
quasi-bound states localized within the cavity. Most 
of the resonances visible in Fig. 0M are also present 
in Fig. [TJa) as nearly horizontal lines - independent 
of the distance between the tips - which suggests that 
they correspond to the quasi-localized stated of the 
cavity. The nearly horizontal lines are very thin in¬ 
dicating a long lifetime of the resonances. The lines 
with a steeper dependence on d t i P in Fig. [7|a) corre¬ 
spond to resonances supported by the tips. The lines 
are wide - indicating a stronger coupling to the chan¬ 
nel. A study probing the conductance as a function 
of the distance between the tips would be an experi¬ 
mental implementation of the stabilization methods! 17 
for detection of localized resonances. 

The zoom of region of the Fig. J7| marked by the 
black rectangle is depicted in Fig. [8fa). In this pic¬ 
ture we can distinguish three types of resonance lines: 
a) vertical - almost independent of the width W of 
the cavity, b) resonances which vary with W and c) 
well visible horizontal resonances - independent of the 
distance between DT which correspond to the quasi 
localized state in the cavity. We found in general that 
the difference between the two first types of resonances 
come from the symmetry of the resonant scattering 
densities. The scattering densities for lines of type a) 
and b) are symmetric with respect to the center of 
the cavity in the x direction. The resonances that are 
independent of W have a nodal surface at the sym¬ 
metry axis, while the other have a maximum on the 


a) 1 G[G 0 ] 0 b) 



Figure 7: a) The conductance of the system depicted in 
Fig. HI as a function of d t i P and cavity width W. Dashed 
rectangle corresponds to the zoomed area shown in Fig. 
[8ja). Red line corresponds to scan with third tip discussed 
in the next paragraph, (b) Conductance of the system 
presented in Fig. [6] as a function of W but without DT. 

symmetry line. The former are strongly localized in 
the center of cavity (see Figs. |8](b-e)) and the latter 
form resonances which are delocalized over the entire 
cavity (see Figs. [8jf-g,i)). The delocalized resonances 
react stronger to the value of W in Fig. 0a). Fig. 
UDM shows a resonance of a third type - of the energy 
that is weakly dependent of du p - which corresponds 
to the resonances supported by the cavity itself. 

In experiments variation of the geometry of the cav¬ 
ity can be accomplished by tuning an electrostatic 
confinement in the y direction due to external gat¬ 
ing. Since the electrostatic confinement potential is 
usually parabolic at its origin we replaced the hard 
wall confinement by 

^parabolic = \u ,2 y 2 = ^ (3meV - u}) 2 y 2 . (5) 

With this potential the effective width of the cavity 
is controlled by the harmonic oscillator energy. The 
parametrization was chosen in order to get better com¬ 
parison with Fig. [7] The value to = 3meV => to' = 
OmeV corresponds to a small value of W in Fig. |b]and 
channel-like behavior of conductance. Value ui = 0 
meV=> uj' = 3meV corresponds to the large values of 
W. Comparing the images Fig. [7] and Fig. [9] we can 
see that both share the same features: a) X-shaped 
resonance lines are present, b) some resonances de¬ 
pend weakly on W as a function of the cavity width 
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Figure 8: (a) Zoom of figure j7^a) marked by the dashed 
rectangle, (b-i) Probability density for the the electron 
incoming from the left for d t i P and W pointed by the ar¬ 
rows. (b-e) Antisymmetric resonances (f-i) symmetric res¬ 
onances. 


(antisymmetric ones) and some vary with du p (sym¬ 
metric states), c) horizontal resonance lines are also 
present. Note that in this case we do not find the 
vertical resonances in the conductance image. This is 
because the variation of ui changes the potential profile 
in whole cavity and not only the width of the cavity 
(as it was in the previous example). 

Concluding this Section, we found that using double 
tip system with the cavity of variable width one may 
tune the system to the specific resonance i) symmetric 
or antisymmetric resonances supported between the 
tips ii) or resonances quasi-localized within the cavity 
(horizontal lines in Figs. [7] and [9]). 


VI. MAPPING THE LOCAL DENSITY AT 
FERMI LEVEL 

Let us now discuss the possibility of mapping the 
LDOS inside the device presented in Fig. [b] In the 
two-dimensional systems the conductance maps can 
be correlated to the local density of states for weak tip 
perturbation when the electron scattering wave func¬ 
tions incoming from the left and right leads are the 
same up to a constant phase. The lat ter ca se can be 
supported^ by the perturbation theorjff^EHI f or w hich 



1 G [G 0 ] 0 
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Figure 9: The conductance of the system depicted in Fig. 
H as a function of dtip and angular frequency lj of the 
oscillator. The parametrization u' = 3meV —co is applied 
in order to keep the same interpretation of y axis (width 
of the cavity) as in Fig. [7|a). 


the first-order correction to the conductance is simply 
proportional to the LDOS. The scattering densities for 
the electron incoming from the left lead plotted in Fig. 
[8])b-i) are highly symmetric since they correspond to 
the resonances induced by the double tip system in a 
symmetric cavity. Thus for inverted current direction 
the scattering density inside the cavity stays the same 
and this should lead good correlation between LDOS 
and corresponding G map. In general we fincPHI that 
the correlation LDOS-G is well resolved near the Fano 
resonances. 

In the following calculation we assume that the dis¬ 
tance between the two tips placed within the channel 
is set to a resonance, and the third tip scans the sur¬ 
face of the cavity for detection of the resonant LDOS. 
The assumed third tip parameters are U t ; p = ImeV 
and Wtip = lOnrn. We need a pronounced variation of 
G as a function of the tip position and the adopted 
height of the tip potential is quite large as compared 
to the Fermi energy kept at 2 meV. 

In Fig. [Iota) we show the conductance of the sys- 
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tem depicted of Fig. [6] as a function of distance be¬ 
tween the two tips inside the channels du p (see Fig. 
[7j) and position of the third tip x t i P along the red line 
in Fig. § We observe a number of resonances for 
a sequence of the tip positions. The resonances ap¬ 
pear in pairs, which have LDOS with a minimum or 
a maximum at the symmetry axis of the device [Fig. 
|10[ c)[ with either a maximum or a zero of LDOS at 
the symmetry axis Xu p = 800 rnn. The du p distance 
which corresponds to resonances can be deduced from 
the conductance dependence of the cavity - double tip 
system (without the third tip) which is plotted in Fig. 
\E d )- In Fig. [l0|b) we present the deviation AG 
from the conductance G' of the system without the 
third tip. Note that AG changes sign in the vicinity 
of the resonance induced by double tip (see the blue 
curve in Fig. [To))d)) , which will lead to a negative cor¬ 
relation of between the LDOS and G map 18 . In order 
to compare the conductance obtained from scan with 
third tip G scan = G(£tip, dtip) with the LDOS on that 
line L scan = LDOS(cc t i P , d t i P ) (for given value of d t i p ) 
we calculate the Pearson correlation r = r(d t i p ) 


T( G scan -<G scan »(L scan (Lscan ))dXtip . . 

A L cr(G scan )cr(L SCQri ) 


where (a) = f a(xti p )dxti p is the average value of 
a, cr 2 (a) = f (a(£ti P ) — (a)) 2 dx ti p is standard devi¬ 
ation of a and Al = 800nm is the length of the scan 
along red line in Fig. [6] Before the calculation of 
r both functions G scan and L scan were normalized to 
range from 0 to 1. The absolute value of the correla¬ 
tion r is plotted in Fig. [Tojd). Note that r increases 
to about 0.9 near the induced resonances which means 
that obtained G map well resolve the LDOS inside the 
cavity. The dips in r{du p ) dependence result from the 
fact that the correlation coefficient changes sign at the 
resonance^]. In Fig. 
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we present the SGM images 
and corresponding LDOS for points A, B and C of Fig. 
|10[ d). Summarizing, the usage of multiple tips allows 
for 1) tuning the system to a resonance and then 2) 
an accurate read-out of the local density of states. 

For evaluation of the LDOS we normalized the scat¬ 
tering densities to the unity within the computational 
box. For off-resonant conditions, where the scatter¬ 
ing density is large only inside the input and output 
channels, we thus obtain a strong reduction of LDOS 
within the cavity. This normalization approach gives 
a close correspondence of LDOS with SGM maps for 
varied du p and Xti P [cf. Fig. 10 . For fixed values of 
these parameters, the Pearson correlation coefficients 


between the maps is not affected by the LDOS nor¬ 
malization. 


VII. TIP AS A PART OF THE SCATTERING 
SYSTEM 


The scanning gate microscope tip can be used as 
a part of the scattering system, with the other tip 
gathering the conductance maps. In particular, the 
tip potential placed near the exit of the quantum 
point contact has been used to form a quantum Hall 
interferometer^. A similar quantum Hall system in¬ 
cluding a potential defect spontaneously formed at the 
exit from the QPC constriction was studied with the 
scanning gate microscopy in RefP 

The SGM with a double tip system can be used for 
observation of the Young interference^ for a set-up 
depicted in Fig. 12 The SGM signal was predicted^ 
to contain signatures of the double-slit interference for 
low-energy transport. The proposed^ set-up requires 
a beam splitter - a central obstacle for the electron 
flow, which may be difficult to introduce in form of 
a fixed gate. A floating one - introduced by the tip 
- is a possible option. In Fig. [12] we plotted a sys¬ 
tem with the QPC slit entering a small open cavity. 
The potential of the tip (schematically shown with the 
gray circle) was introduced near the exit of the cavity 
to an unconfmed half-plane. The scattering density 
current \j(x,y)\ plot in Fig. 12 shows formation of 


the interference pattern characteristic for the Young 
interference^. The other tip used as an electron flow 
detector could then be used for readout of the double¬ 
slit interference in the SGM conductance images. 


VIII. CONCLUSIONS 

We have studied numerically the possible applica¬ 
tions of the multitip scanning gate microscopy on the 
systems with 2DEG. We found that the measurements 
using a pair of tips allows for 1) readout of the Fermi 
wavelength 2) detection of the potential defects with 
indication of their position and profile, 3) filtering of 
specific transverse modes 4) tuning the system into a 
resonance that allows for a reliable detection of the 
local density of states with a third tip. 
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Figure 10: a) Conductance G of the system as a function of third tip position *ti P (moved along the red line in Fig. [6| 
and distance between the two tips (red line in Fig. [7](a)) b) The same but AG — G — G' is plotted, where G' is the 
conductance of the system without the third tip c) is the LDOS. d) (blue line) Conductance G' of system as a function 
of dtip without third tip. (red line) Absolute value of Pearson correlation r between (a) and (c) images calculated from 

Eq. 0 



x Up [nm] 


Xu p [nm] 


[nm] 


of Krakow Smoluchowski Scientific Consortium from 
the funding for National Leading Reserch Centre by 
Ministry of Science and Higher Education (Poland). 


Figure 11: a-c) The SGM image d-f) and LDOS obtained 
for point A,B and C in Fig. |10[ d). Point A corresponds to 
the symmetric resonance induced by the double tip system 
in work point where the correlation r is negative. Point B 
corresponds to the antisymmetric resonance with positive 
correlation r. Point c corresponds to a zero of T in the 
absence of the tip. 
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